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The transverse polarization of forward A hyperons produced in high-energy p-A collisions is ex- 
pected to display an extremum at a transverse momentum around the saturation scale. This was 
first observed within the context of the McLerran-Venugopalan model which has an i-independent 
saturation scale. The extremum arises due to the fct-odd nature of the polarization dependent frag- 
mentation function, which probes approximately the derivative of the dipole scattering amplitude. 
The amplitude changes most strongly around the saturation scale, resulting in a peak in the polar- 
ization. We find that the observation also extends to the more realistic case in which the saturation 
scale Qs is x dependent. Since a range of x and therefore Qs values is probed at a given transverse 
momentum and rapidity, this result is a priori not expected. Moreover, the measurement of A 
polarization over a range oi xp values actually provides a direct probe of the x dependence of the 
saturation scale. This novel feature is demonstrated for typical LHC kinematics and for several phe- 
nomenological models of the dipole scattering amplitude. We show that although the measurement 
^ , will be challenging, it may be feasible at LHC. The situation at RHIC is not favorable, because the 

O ' peak will likely be at too low transverse momentum of the A to be a trustworthy measure of the 

saturation scale. 
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I. INTRODUCTION 



It is well-known that A hyperons produced in collisions of unpolarizcd hadrons are to a large degree polarized 
perpendicularly to the production plane. Even though the origin of this phenomenon has not been clarified 
fully yet, for sufficiently large transverse momentum pt of the A, one expects that a parton description must be 
\ applicable. In Ref. [l| it was shown that the available A polarization data for pt > 1 GeV can be described within 
a factorized approach by employing a polarization and transverse momentum dependent fragmentation function. 



' This function, denoted by D, describes the fragmentation of an unpolarized quark into a transversely polarized 



A and is an odd function of the transverse momentum kt of the quark w.r.t. the A momentum This kf-odd 
nature implies that it is essentially accompanied by the first derivative of the partonic cross section w.r.t. kt, 
unlike the unpolarized A fragmentation function D, which is kf-even. This turns the observed A polarization into 
a useful tool, which in this paper will be applied to the study of the x-dependence of the saturation scale. 
I Once the polarization dependent fragmentation function is known from data in the large x or DGLAP region, 

polarization measurements in other kinematic regions could point to changes in the underlying physics. In Ref. 
[3| this was discussed specifically for the saturation region, the region of small momentum fraction x where the 
gluon density is very high and is expected to saturate according to the nonlinear evolution equations of relevance 
^ ' in that region. The saturation region is characterized by the so-called saturation scale Qs, i.e. the momentum 
5^ , scale at which saturation effects become sizable. It was noted within the context of the McLerran-Venugopalan 
model which has an x-independent saturation scale Qs, that the negative valued A polarization displays a 
minimum at a transverse momentum approximately equal to Qs- Here, we want to investigate the more realistic 
case in which the saturation scale is x dependent [1,[1,0]. It may be expected that, since now a range of Qs-values 
is probed, the minimum of the polarization is smeared out, and possibly not recognizable anymore. However, in 
this paper we will demonstrate that this is not the case. In fact, the pronounced minimum of the polarization can 
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even be used to probe the x dependence of the saturation scale. This makes the observable of potential interest 
for collider experiments at RHIC, LHC and a future electron-ion collider, the EIC. 

Having discussed this promising use of A polarization, let us now address the possibilities of measuring it 
in the small- a: region. In high energy scattering the polarization of a spin- 1/2 final state hadron can usually 
only be measured through self-analyzing parity violating decays. Exploiting this property it was demonstrated 
already more than 30 years ago [1, 0| that A hyperons emerging from unpolarized p -A collisions are polarized 
perpendicularly to the production plane (cf. Ref. [13] for an extensive review of data). In the fixed target 
experiments performed at typical center of mass energies ^/s 20 GeV, the transverse momentum dependence 
of the degree of polarization shows the characteristic feature that after a linear rise up to Pi ~ 1 GeV, it stays 
remarkably constant up to the highest measured values pt ^ 4: GeV. This behavior was found to be independent 
of the specific values of and atomic number A. For larger pt values one expects the asymmetry to fall off as 
1/pt, but this has not been observed yet. None of the measurements performed thus far are in a kinematic region 
where the target could be considered dense, i.e. in the saturation region. In Ref. [1] it was pointed out that the 
characteristic flat behavior observed for pt values of a few GeV, would in that case no longer be present, but rather 
an extremum should be visible, located at pt ~ Qs. In other words, the observed plateau should turn into a peak 
as saturation effects set in and Qs becomes a relevant scale. Since Qs grows with 1/x and A, one expects this to 
happen when ^/s and/or A are increased significantly. In addition, it also helps to consider large rapidities, i.e. 
forward A production, in order to decrease the x values probed. At RHIC forward A's with rapidities of around 
4 would begin to probe the small- a; region according to a dipole scattering description [Tl!|. The possibilities to 
probe small x-values are of course greater at LHC, where due to the much higher ^/s much less forward A's are 
required. For completeness, we recall that at RHIC d-Au collisions have been studied at energies of ^ 200 GeV/ A 
in the nucleon-nucleon center of mass frame. At LHC p-Pb collisions will be performed at ^^snn — 8.8 TeV, but 
these do not take place in the nucleon-nucleon center of mass frame which leads to a rapidity shift from lab frame 
to center of rapidity frame of about half a unit. In principle also the p-p collisions at LHC are of interest here, 
due to the very large energy: ^/s = 14 TeV. 

Experimentally the measurement of forward A's and their polarization may be hampered by the often restricted 
particle identification capabilities in the forward region. Two-thirds of the time A's decay into protons and 
negatively charged pious: A -~^ pi:^ . The angular distribution of the decay in the A rest frame is used to 
determine the polarization of the A. Unfortunately, protons are usually hard to identify in the forward region. In 
that case the only alternative may be to use that the A's decay one third of the time into neutrons and neutral 
pious (and subsequently, two photons): A — > titt^ ^ 7177. Neutrons, tt^'s and photons have been identified in the 
forward region at RHIC, hence this alternative may be feasible [l2| and may in fact be the only way of measuring 
A polarization in the forward region at RHIC, LHC or EIC. We will proceed with our investigation under the 
assumption that forward A polarization measurements will be possible in the future. 

The outline of this paper is as follows. In section [TT] we discuss the A polarization asymmetry in terms of 
the relevant polarization dependent fragmentation function and the dipole scattering amplitude. This discussion 
repeats the essentials from Ref. Q in order to set the notation and to explain why an extremum is expected at 
Pt Qs, but also it includes the details of various phenomenological models for the dipole scattering amplitude 
that were considered in the literature. As mentioned, in Ref. ^] only the McLerran-Venugopalan model was 
considered, but here we will focus on more recent models that employ an x-dependent saturation scale. In section 
mil we discuss model results for the A polarization observable for LHC kinematics mainly and point out the 
generic qualitative features. Achieving realistic quantitative predictions for the degree of A polarization will not 
be our aim, due to the large uncertainty in the polarization dependent fragmentation function. Nevertheless, an 
estimate can be given of the range oi xp values required to observe the x dependence of the saturation scale, as 
the Pt dependence of the A polarization is found to be less model dependent than its absolute value. Prospects 
for RHIC are also briefiy discussed, but no results will be shown. The reason for this is that at RHIC the peak 
is likely situated at a pt below 1 GeV, where the considered framework would not be appropriate. We end with 
conclusions. 
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II. TRANSVERSE A POLARIZATION DESCRIPTION AT SMALL x 



According to Ref. the transverse polarization of forward A's produced in unpolarized -p-A collisions is 
approximately given by 



Jxf 



dx X 



E,/,/p(x,Mpi^A/, [^puX,)+f,/,{x,f,})DA^, {^Pu^^ 

(1) 

where yh and pt are the rapidity and the transverse momentum of the produced A, xp — PtI \ls exp[2//i], and x 
and xi = cj^jix s) are respectively the momentum fractions of the parton in the proton and the heavy ion (referred 
to as the "target"). Note that three different values of xi enter in this polarization expression, in conjunction 
with the three different transverse momenta qt = x/xpipt ± ^) and qt = x/xpPt of the scattered parton. The 
expression ([T]) is based on the asymmetry expression of Ref. [l| combined with the dipole picture description of 
the cross section of Ref. [ll|. In the dipole formalism, Np describes the scattering of a quark off a nucleus, while 
Na describes gluon-nucleus scattering. For details of the derivation and justification of the approximations we 
refer to Ref. j3|, where the dipole scattering amplitude Np is denoted by C. 

The polarization dependent fragmentation function A^D^t /q is parameterized in terms of the unpolarized one 



Da/^ of Ref 



where 



DKyq{z,n^) = f^{z)Dp,/^{z,n^) 



^ ^Ngz^^il- zf^ , Nd = -28.13, TV, = 57.53 , c, = 11.64, = 1.23, 

and the average transverse momentum is given by 

k^{z) = 0.66^"-^^(l - z)°-5 GeV . 



(2) 



(3) 



(4) 



We emphasize that there is a large uncertainty in this parameterization extracted from fixed target data [l| , so that 
the numerical results presented below should only be viewed as qualitative, not as quantitative predictions. Future 
collider data from LHC could be used to obtain a more trustworthy parameterization, for instance through the 
A+jet observable recently pointed out in Ref. [14], which deals with A's at midrapidity where particle identification 
does not pose a problem. 

Using the McLerran-Venugopalan (MV) model for the dipole scattering amplitude it was shown in Q that 
the Pt distribution of the transverse polarization displays a peak that is directly related to the saturation scale 
Qs- However, the MV model does not incorporate evolution in x; Qs is constant. Here we want to investigate 
whether Va as a function of pt still possesses an observable peak when described using a more realistic dipole 
scattering amplitude including x-evolution. A phenomenologically successful model with an x-dependent Qs is 
the Golec-Biernat and Wiisthoff (GBW) model [HI, which was able to describe small-x DIS data^. Here we will 
focus on two different modifications of the GBW model that were introduced to describe RHIC d-Au data: the 
DHJ model and the geometric scaling (GS) model of Ref. ^1^. Both are of the form 



NF{qt,X2) = / (fn e 



1 — exp — 



4 2 
9 * 



(5) 



where Na is obtained from Np by replacing irf (5^(0:2) by (5^(0:2), and the saturation scale is given by [iJ 



,(X2) = IGeV 



(s) 



A/2 



(6) 



^ Note that the GBW model was found to be inconsistent with newer, more accurate data and requires some modification at larger 
Q^, see for example Ref. ,161] . 
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where the parameters xq ~ 3 x 10"'' and A ~ 0.3 were fitted to the small- a; DIS data. In the GBW model the 
so-called anomalous dimension 7 is simply equal to 1. The DHJ model incorporates expectations on the behavior 
of 7 from BFKL/BK evolution [13, namely a logarithmic dependence on qf /Q'^{x), and violations of geometric 
scaling that are proportional to 1 /y (at large y) 

. (n .\ ^ |l0gfaVQK^2))| 

+ dy^ + \ ^og{qt /Qi{x2))\ 

where 7^ — 0.6275, y = logl/a;2 is minus the rapidity of the target parton. The saturation scale Qs{x) and the 
parameter A are taken from the GBW model, as given in Eq. and d = 1.2 was fitted to data. Here Qs includes 
the additional factor A^^^, where for large atomic numbers A usually a lower, effective number A^e is used to 
account for impact parameter dependence. 

The parameterization of 7dhj, which is based on the one given in Ref. '20'|, is well motivated by expectations 
from small-x evolution that are valid only for qt > Qs @, [2l|. In contrast, the continuation to the saturation 
region qt < Qs is rather undetermined. In the case of hadron production at RHIC [13, [lB| this is not crucial since 
this region is hardly probed. But the polarization observable discussed here is sensitive to 7 around Qs, not only 
to qt > Qs- Hence, the continuation of 7 to the saturation region affects the polarization around the peak. 

The DHJ model was found to describe well forward hadron production in d-Au collisions at RHIC, but it fails 
at central rapidity [l^. In particular, the logarithmic rise of 7 proved to be too slow to describe the larger x 
central rapidity data. Also the scaling violations of the DHJ model could not be resolved in the data. In order to 
investigate to what extent the RHIC data establish the small- a; properties incorporated in the DHJ model, a new 
model was put forward that is not only exactly geometrically scaling, but in addition features a stronger rise of 
7 M 

7GsH=7s + (l-7.) (J,_^^;^ . (8) 

Here, a = 2.82 and h = 168 were fitted to the data. This GS model turned out to be able to describe the d-Au 
data well at all rapidities. It must be emphasized though that when restricted to the forward data, or equivalently 
the smaller x data, both models describe the data equally well. From the comparison of the model predictions to 
future LHC p-Pb and p-p data one should be able to learn which rise is more appropriate at small x. Given this 
uncertainty, here we will use both models to study the transverse A polarization. 

To shed light on the peak in the pt distribution, we will separate Va into a pt-dependent and an x ^-dependent 
part in the following way. To good approximation the integrals in ^ are dominated by a value of xp/x = z 
that is independent of pt and only moderately dependent on a;^. Due to the large power Cq, which suppresses 
small ratios z in the numerator, the values of z effectively probed in the numerator and the denominator are 
different. We will denote the value that dominates the numerator with z and the smaller one that dominates the 
denominator with z' . Of course in the kinematic limit xp ^ \, both z and z' must become equal to 1. In the 
following analysis we will stay away from this limit and assume that xp stays smaller than roughly 0.5. Ignoring 
the gluonic contributions, which is a good approximation when xp is not too small, we can approximate ([1]) in 
the following way 



Vk^Pi^xp) 



Eg -Pa/9 (^) {xfIz) fq/p{xF/z, f^{z) 

D^/qiz') {Xp/z') fq/p{xFlz', fi'j) 



Xp z' s 



(9) 



Since z and z' are considered constant, Eq. Q now depends on pt through the function Np only. This is true 
assuming the factorization scale /z/ to be constant. Below we will mostly choose Hf = pt though, but this will 
turn out not to make much difference. We further note that since is only around 0.3 GeV or smaller for all 
relevant values of z, we can expand Np {\{pt — ^t')) — Np (^(pt + k^ )^ in terms of k^/pt, requiring pt > I GeV 
throughout this paper, 

Np - k?)) - Np {^{pt + fc?)) « -2 ^ ^ . (10) 

z dqt 
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Here we have suppressed the exphcit dependence on X2 for convenience and we will do so frequently below. Writing 
the dipole scattering amplitude in terms of a dimensionless function Np, 

NF{qt,X2) = —Nf{w = qt/Qs{x2),X2) , (11) 

we can express Eq. (fTU|) in the following way, 

1.0 o_ 

Nf (lipt - fc?)) - Nf + fc°)) « 2^^{2Nf{w) - wN'piw)) . (12) 

Using this result, we can split off the pt-dependence of the transverse polarization and write 

VAiPt Xf) ~ -PA/gl^) {xf/z) fg/p{xF/z,fl}) f^{z)/z fcO ^^^^ ^^^^ 

J2qDA/q{z'){xF/z') fq/p{xF/z',flj) Q s z' 

where we have defined the pt-dependent part of Va as a separate function F{w,w'), 

F(^,^0^1^^^(-)""^^("). (14) 
w NFiw') 

From the asymptotic behavior of F it can be seen that it must have an extremum. From Eqs. ^ and ([TT|) . it 
follows that Nf oc for large w, and hence that F{'w,w') will approach 2(1 + 7)/^. On the other hand, in 

the deep saturation regime the function pT|) is proportional to w"^, so that F{w, w') vanishes as w — > 0. Therefore, 
the function F{w,w'), and consequently Pa, must have a peak as it connects these two asymptotic behaviors. 
Without saturation there could also be a peak in Pa, but one would in that case not expect the extremum to be 
rather sharply peaked at a perturbative scale of a few GeV. Such a peak would be a sign of saturation, especially 
if it increases towards larger transverse momenta with increasing energy. Of course, there could be a plateau-like 
extremum, as it appears to be the case at low energies. However, the MV model calculation of Ref. Q clearly 
shows there to be a pronounced peak, with a position proportional to the (constant) saturation scale Qs- If this 
proportionality holds when Qs evolves with x, the location of the peak in pt would be a direct probe of the running 
of Qs through its dependence on xf- If however the peak position depends also explicitly on xf, the running of 
Qs cannot be reconstructed from the peak position. Because the probed values of z' /z — w/w' depend on xf, 
this means that we have to check that they do not influence the position of the peak. Fig. [1] shows F{w,w') 
for various values of values of w/w' = z' / z ranging from 0.25 to 1, using a dipole scattering amplitude with a 
constant 7 = 0.6275. The curves indeed have a clear maximum^ as a function of w. The position of the peak 
hardly depends on w'/w if w' /w is not too close to 1, i.e. away from the kinematic limit — > 1- Hence, we 
conclude that the peak of F is located at an approximately constant value of w. As mentioned, this means that 
the minimum of Pa docs not explicitly depend on XF^ but only through the saturation scale Qs- 

We find that all this remains true not only for different constant 7's, but also for the DHJ and GS models. 
The GS model actually leads to the same peak position as constant 7 = 7^, because it differs only little from 
lGs{w = 1) = 7s in the saturation region w < 1, where the peak of F is located. The DHJ model gives a slightly 
different peak, one that depends on the continuation into the saturation region qt < Qs, as will be discussed 
further below. 

One can estimate the xf dependence of the peak of the resulting pt-distribution as follows. Since the peak of F 
is at a constant value ofw = zpt/Qs, where z is roughly constant as well, the peak in pt is directly proportional to 
Qs{x). Because the dominant value of x that is probed depends on both pt and xf, the peak position Pt°^^ will 
depend on xf- As the probed value of z = xf/x in the integrals in Eq. ([1]) is to good approximation constant, 
the target momentum fraction X2, which sets the saturation scale Qs{x2), is given by 

X Tt^ 1 Tt^ 

X2 = X exp{-2yh) = ^ — OC -. (15) 

xi, s Xf s 



^ Since the polarized part of the fragmentation function is negative for the u and d quarks that lead to the dominating contribu- 
tions, "Pa will have a negative valued minimum, which for convenience will sometimes also simply be referred to as a peak. 
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FIG. 1: The function F{'w,w') for 7 = 0.6275 and various ratios w/w' 



Using this relation, we can estimate the a;F-dependence of the peak position p^'^^^ of Va- Assuming that the 
saturation scale is given by a power law in 1/x, Eq. ([6|), we see that 

\ A/2 

pP"" oc QsixF^p^""^) (X Qo ( (16) 




=^ pr\xF) cx Qo4'/' (^) ^'=lf\ ^^^^ 

Hence, we conclude that the running of the peak position with is a clear indication of the running of the 
saturation scale Qs(a^2)- Moreover, the power A can be reconstructed from the behavior of the peak position as a 
function of xp- 

The parameterization of Qs in Eq. © is not just based on the GBW model fit to DIS data. The specific power 
law dependence on 1/x stems from theoretical arguments. Small-x evolution equations, such as the GLR [23 |. 
BFKL [13 and BK [23| equations in the fixed coupling constant case, result in such a behavior and determine 
the power A (typically they yield A « 0.9). In the running coupling case the functional form of Qs is different. 
But over the limited range of experimentally accessible values of x, Qs can still be approximated by a power 
law like behavior. In this way the specific value of A = 0.3 that best describes the DIS data, can be effectively 
accomodated. This implies that A may be different in p-A collision where a different kinematic range is probed. 

In the discussion above we have expanded Np in terms of k^/pt requiring pt > 1 GeV, so that pt is in the 
perturbative regime. This should be considered a minimal requirement for the present dipole description to be 
applicable and also from the perspective of the scale choice = pt it is a sensible lower bound. Therefore, below 
we will only discuss results for which p^"^^ is in the perturbative regime. 



III. TRANSVERSE A POLARIZATION RESULTS 



Here we will present our numerical estimates of the transverse polarization ([T]). For the fragmentation functions 
we have chosen the leading order (LO) functions given in [l^ and for the parton distributions the CTEQ5 LO 
ones [2^. In [3] it was shown that certain effects of higher order can be taken into account by the DGLAP 
evolution of fragmentation functions and parton distributions at the scale set by pt- Unless stated otherwise, we 
will therefore set the factorization scale to /i/ — pt- We will return to the fXf dependence of the results later on. 
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We first discuss the pt distribution of Va for constant values of xp- Here we will first give the result for p-Pb 
collisions at LHC at y/s = 8.8 TeV explicitly and later point out how they compare to p-p collisions at LHC and 
d-Au collisions at RHIC. For the saturation scale we will use the GBW parameterization ^ with Qq = 2.7 GeV^ 
instead of 1 GeV^ by taking Agg = 20. Figure[2]shows the resulting Va- It has been calculated for dipole scattering 
amplitudes with various constant values of 7 from 0.5 to 0.9. The increasing magnitude of the polarization with 
increasing is due to the polarized part of the fragmentation (for larger x p larger values of z are probed) . 
The anticipated rise of the peak position with xp can be clearly observed. Furthermore, the peak position rises 
approximately linearly with 7 for all considered values of xp. 



XF = 0.1, 7 = 0.5-0.9 

XF = 0.3, 7 = 0.5 - 0.9 




FIG. 2: Vh for various constant 7 for p-Pb collisions at ^/s = 8.8 TeV. Curves for smaller 7 have their minimum at smaller 
Pt- 

Figure [3] shows the polarization for various a;;' as a function oi pt, but now for three 7's that are all equal at 
the saturation scale: a constant 7s, 7dhj and 7gs- As expected, the difference between the polarization for 7s 
and 7gs is very small because in the saturation region 7gs differs only mildly from the value 703(9* = Qs) = Js- 
Because of the uncertainty in the continuation of the DHJ parameterization in the saturation region, the estimate 
of Va around the peak is ambiguous. If we would for instance continue 7dhj by keeping it constant for qt < Qs, 
we would obtain roughly the same result as for 7gs- The fact that the DHJ and GS models yield similar results 
for the observed behavior of the peak indicates our findings are rather robust and to a certain extent model 
independent. In contrast, the magnitude of the polarization is subject to considerable uncertainty, mostly due to 
the parameterization of D, but also somewhat due to the choice of the factorization scale. 

Thus far we have used a factorization scale fif — pt- However, in Ref. Q — Qs was considered, which may 
also be a natural choice. In the present case that would lead to an x-dependent factorization scale. Figure [H 
shows Va with 7gs for three different factorization scales, ^/ = pt, jJf'f = Qs, and a constant scale /i/ = 1 GeV. 
As can be seen, for constant x p the shape of the pt distributions is rather independent of the factorization scale. 
The normahzation does depend on the choice oi Hf, but still only moderately. Therefore, choosing fif = pt does 
not noticeably affect our claim that the x p dependence of the peak momentum directly probes the x dependence 
of the saturation scale. 

Fig. [5] shows the xp dependence of the peak position of the pt distribution for various choices of 7. The lines 
for constant 7 confirm that the peak position scales linearly with 7. Moreover, for not too large xp, the power 
law rise of Pf°^^ with xp is consistent with the result we obtained in Eq. (|17p . including the fact that the power 
is independent of 7. As expected, the results for 7gs(w) (0 and the constant 7 = 7s are very close to each other. 
The curve for 7dhj is similar to that of a constant 7 that is slightly larger than 7s. This is because 7dhj rises 
rather quickly in the saturation region as qt decreases. 
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FIG. 3: Va in p-Pb collisions at ^/s = 8.8 TeV, for 7dhj, 7gs and a constant 7^. The top lines correspond to xf — 0.1, 
the lowest to xf — 0.5. 
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FIG. 4: Va in p-Pb collisions at LHC for the scaling 7gs (O and three different choices of the factorization scale fif — pt,Qs 
and 1 GeV. The top lines correspond to xf — 0.1, the lowest to xp ~ 0.5. 



As can be seen from Fig. [5] all slopes are also numerically consistent with the power A'/2 in Eq. lfT7|) for A = 0.3. 
This implies that an increase in xp by a factor 5 leads to a shift in the peak position of approximately 20%. 
This can be seen directly in Fig. |3]too, when comparing the peak position at xp = 0.1 and 0.5. Unfortunately, 
this is not a large shift, but it does give an estimate for the precision with which the peak position needs to be 
determined. It should be mentioned that, as discussed before, the value of A may be different in p-A collisions 
than the one taken from the analysis of the small- a; DIS data. A larger value of A would of course result in a 
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stronger x dependence of Qs and therefore in a larger xp dependence of the peak position that would be easier 
to observe. At small xp, where the position of the peak is less pronounced, it will be harder to determine than 
at large xp. The value of pt/Qs at which the peak is situated depends -too good approximation linearly- on 7 in 
the saturation region qt < Qs- As can be seen in Fig. [31 the peak is located for 7gs and 7 = 7^ at almost the same 
postion. We find empirically that in these cases the minimum in the pt distribution shows up at w 0.55, i.e. 
Pt ~ 0.55 2: Qs, where z rises slightly with xp from 0.7 to 1 in the limit xp —> 1. Depending on the continuation 
of the DHJ model to the saturation region the peak is situated at a different w. For the continuation ([7]) the peak 
shows up at Pt ~ 0.60 2: Qs, since it rises again towards smaller qt in the saturation region. 

Similar results for Va are obtained for p-p collisions at LHC and d-Au collisions at RHIC. Again an extremum is 
observed at around one half times the saturation scale that shows the same xp dependence as for p -Pb scattering 
at LHC. However, due to the different kinematics and targets, Qs and hence the position of the peak Pt^'^^ is 
in both cases lower. Following the same line of arguments leading to the xp dependence of Pt^'^^ (|17p . one can 
estimate its ^/s and Qq dependence, 

pr\xp, V7, Q'o) = pr\xp, v~s, Qo) ■ (18) 

\QoVs / 

For p-p at LHC Qo is 1 GeV and ^/s = 14 TeV. Hence, the peak position is expected to be reduced by a factor 
of 1.3 with respect to p-Pb collisions at i/s = 8.8 TeV. An explicit calculation confirms that this estimate works 
very well, i.e. the x^-dependent extremum is expected to show up approximately between 1.5 and 2.0 GeV for 
Xp ~ 0.1 — 0.5. For d-Au collisions at RHIC the probed values of X2 are less small due to the smaller energy. 
Hence, the probed values of Qs and hence pP°^^ are reduced even more, namely by a factor of 2.4 compared with 
p-Pb collisions at LHC, which may situate it below the perturbative regime pt ^ 1 GeV, even for constant values 
of Xp = 0.1 — 0.5. However, given the uncertainties in e.g. the values of Qo and A, a peak in the perturbative 
region is not ruled out, especially for larger xp. From this perspective it may still be worthwhile to investigate 
this observable at RHIC. 

Up to now we focused on the calculation of Va at constant xp, where the dependence on is not that large. 
However, from an experimental point of view it might be more convenient to measure Va for constant rapidities 
Uh- As demonstrated before, there is a clear peak in the pt distribution at fixed xp, which is at different locations 
for different xp. Therefore, since at fixed yh a range of values of contributes, the peak will be smeared out 
to some extent (this can also be observed for the DHJ model predictions of single spin asymmetries in forward 



10 



pion production in the collisions of transversely polarized protons with unpolarized protons [26l|). Hence, it is not 
clear a priori whether the peak remains observable and whether the peak position is still a clear probe of the 
saturation scale. 

For LHC kinematics, we know from the previous analysis that a peak at transverse momenta larger than 1 GeV 
requires xp = Pt/\/s exp[?/^] ^ 0.01. At LHC such a peak is thus only expected in the forward region yh 4. 
Fig. [6] shows for p-Pb scattering at LHC at ^Jsnn — 8.8 TeV, for values of yt = 4, 5, 6. Indeed, the extremum 
is in these cases located at a pt larger than 1 GeV, but it is much less pronounced than at fixed xp and for 
the GS model less recognizable than for the DHJ model. We also note that the magnitude of the asymmetry is 
considerably reduced compared to the fixed xp case. 



-0.005 



-0.01 
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FIG. 6: Va in p-Pb collisions at -^s = 8.8 TeV, for constant using 7gs and 7dhj- The top lines correspond to yh ~ 4, 
the lowest to yn =6. 

At RHIC the saturation scale becomes roughly of the order Qs>l GeV for forward A's with rapidities of around 
4. However, unlike for the MV model, the peak position for the DHJ and GS models is located considerably below 
Qs, that is, below = 1 GeV. An explicit calculation of Pa for RHIC confirms that even for yh — i a, peak is not 
expected to be above 1 GeV. In other words, Qs{x) can presumably not be extracted in a trustworthy manner 
from a fixed y^ = 4 analysis at RHIC, unless Qq and/or A turn out to be larger than expected at present. 



IV. CONCLUSIONS 



The transverse polarization of A particles displays a peak at the saturation scale when described using the MV 
model for the dipole scattering amplitude. We find that in the more realistic case where the dipole amplitude 
depends on x, such a peak in the pt distribution remains. The position of the peak, Pf°^^, is still proportional 
to Qs, and therefore offers a direct experimental probe of this scale. For fixed values of xp, the x dependence of 
Qs can be reconstructed from the xp dependence of p^'^^^. It would be very interesting to compare the function 
Qs{x) obtained in p-A collisions in this way with the GBW model one that was obtained from DIS data, in order 
to establish consistency among the descriptions of all available data. The power A in Qs ~ determines how 

strongly the peak varies with xp. Using A — 0.3 as obtained from DIS, which according to a dipole scattering 
description is compatible with forward hadron production d-Au data of RHIC, we have obtained the following 
results. In p-Pb collisions at LHC, for values of xp that are between 0.1 and 0.5, the position of the peak is 
expected between pt = 1.5 and 2.5 GeV. This result is obtained for a range of dipole models that includes the 
DHJ and GS models. In p-p collisions, the position of the peak is reduced by a factor of 1.3, but is still in the 
perturbative regime. In d-Au collisions at RHIC, the position of the peak is smaller by a factor of 2.4 with respect 
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to p-Pb at LHC, due to the much smaller energy. Hence, observing the peak in the perturbative regime at RHIC 
seems unlikely, except perhaps at even larger xp values. 

For fixed values of the rapidity instead of xp, the peak in the pt distribution gets smeared out and is reduced 
in size. Moreover, in this case the polarization peaks in the perturbative regime ^ 1 GeV only for A rapidities 
of 4 or larger in p-Pb collisions at LHC. Therefore, A polarization LHC data at fixed xp are best suited for the 
purpose of establishing the x-dependence of Qs in p-A collisions. 

Even though the presented quantitative estimates are to some extent model dependent, the qualitative features 
of the A polarization, i.e. the position of the peak with respect to Qs and its running with xp, are expected to 
be generic for the small-a; region. This offers a unique possibility to probe Qg directly in p-A collisions. 
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